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The mechanism of heat transfer and the contribution of electron-phonon coupling 
to thermal conductance of a metal-semiconductor interface remains unclear in the 
present literature. We report ab initio simulations of a technologically important 
titanium silicide (metal) - silicon (semiconductor) interface to estimate the Schot- 
tky barrier height (SBH), and the strength of electron-phonon and phonon-phonon 
heat transfer across the interface. The electron and phonon dispersion relations of 
TiSi 2 with C49 structure and the TiSi 2 -Si interface are obtained using hrst-principles 
calculations within the density functional theory (DFT) framework. These are used 
to estimate electron-phonon linewidths and the associated Eliashberg function that 
quantihes coupling. We show that the coupling strength of electrons with interfa¬ 
cial phonon modes is of the same order of magnitude as coupling of electrons to 
phonon modes in the bulk metal, and its contribution to electron-phonon interfacial 
conductance is comparable to the harmonic phonon-phonon conductance across the 
interface. 
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I. INTRODUCTION 


Thermal transport considerations form an increasingly important element in the de¬ 
sign of miniatnrized electronic devices. All semicondnctor devices possess metal contacts, 
and hence the stndy of transport throngh metal-semicondnctor interfaces is a technolog¬ 
ically relevant problem. The stndy of fnndamental mechanisms of heat transfer across a 
metal-semicondnctor interface has received significant attention over the last two decades^^-. 
Stoner and Marisi>^ measnred interfacial thermal condnctance between diamond and sev¬ 
eral metals (Ti, An, Al, and Pb) nsing picosecond optical techniqnes. They reported an 
anomalously high thermal interfacial conductance that exceeded the maximum thermal con¬ 
ductance possible from harmonic phonon transmission across the interface. Stevens et ah- 
measured the thermal conductance of a series of metal-dielectric interfaces involving Cr, 
Al, An and Pt on Si, sapphire, GaN and AIN substrates. The experimental measurements 
were compared with predictions of the diffuse mismatch model (DMM)^, which under es¬ 
timated thermal conductance for impedance-mismatched interfaces while it over estimated 
the thermal conductance of well matched interfaces. 

Electrons are the primary heat carriers in a metal while phonons or lattice vibrations are 
the primary carriers in a semiconductor. At the interface between a metal and a semiconduc¬ 
tor, an energy transfer channel from electrons in the metal to phonons in the semiconductor 
must exist. Numerous theoretical models have been developed^^— to predict the interfacial 
thermal conductance of a metal-semicondnctor interface. These models broadly fall into two 
categories as shown schematically in Figure [H In Mechanism A (see Figured]), electrons 
in a metal couple only with phonons in the metal, and these phonons transfer energy to 
phonons in the semiconductor across the interface. The total interfacial thermal resistance 
can be calculated using a series resistance model as explained in Refs. 0 and y. Majumdar 
and Reddy^ used a two-temperature model coupled with DMM to predict total interfacial 
thermal conductance. The volumetric electron-phonon coupling constant G was estimated 
by fitting to experimental data, and the inclusion of electron-phonon coupling resistance 
to the DMM resistance improved the agreement of the theoretical model with experiments. 
More recently, Singh et ah- obtained the electron-phonon coupling constant directly from 
Fermi’s golden rule by assuming a deformation potential for electron-phonon scattering. The 
magnitude of the deformation potential was obtained by htting to the electrical resistivity 
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of bulk metal. 



FIG. 1. Two mechanisms of heat transfer across a metal-semiconductor interface. Mechanism A 
is a two-step energy cascade where electrons in metal transfer energy to phonons in the metal and 
only phonon-phonon energy transfer occurs across the interface. In Mechanism B, electrons in 
metal couple directly with lattice vibrations in the semiconductor and this occurs in parallel with 
the phonon-phonon energy transfer across the interface. 


In Mechanism B (see Figure [T]), a direct coupling exists between electrons in the metal and 
phonons in the semiconductor. Huberman and Overhauser- were among the hrst to propose 
such a mechanism to explain the anomalous high interface conductance reported for a Pb 
- diamond interface^. They used a deformation potential model to estimate the coupling 
between metal electrons and interfacial phonon modes. Sergeev^*^ also modeled inelastic 
electron-phonon scattering at the interface using an analogy with inelastic electron-impurity 
scattering. More recently, Mahan^ derived an analytical theory for predicting thermal 
conductance between a metal and an ionic insulator. The scattering matrix elements were 
derived taking into account the interaction between conduction electrons in metal and the 
potential generated by images of ions in the insulator. 

High-hdelity techniques such as the atomistic Green’s function (AGF) method and molec¬ 
ular dynamics exist for determining the phonon-phonon coupling resistance across an inter¬ 
face. However, models for electron-phonon coupling near a metal-semiconductor interface 
assume a particular mechanism of transport (Mechanism A or B as shown in Figured]) and 
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derive interface conductances using deformation potentials as empirical parameters. Such 
models offer insights into the various heat transfer mechanisms associated with electron- 
phonon coupling, but they do not possess the hdelity to assess the different heat transfer 
pathways shown in Figure 1 quantitatively. A unihed approach to obtain a quantitative 
estimate of the thermal conductances associated with the various energy transfer pathways 
(electron-phonon and phonon-phonon) in Figure [T] is missing in the existing literature. 

In this paper, the electron-phonon coupling strength is determined directly from hrst- 
principles density functional theory (DFT) calculations. We obtain the strength of electron- 
phonon coupling for a bulk metal and for a metal-semiconductor interface using electron- 
phonon coupling linewidths derived from hrst-principles calculations. Relevant limits of 
phonon transport across the interface are treated in a simplihed manner using the radiation 
limit and diffuse mismatch models. The radiation limit gives an upper limit of harmonic 
phonon interface conductance and enables us to assess the contribution of electron-phonon 
coupling relative to the contribution of phonons to the interface conductance. 

Transition metal silicides are widely used as metal contacts in silicon integrated circuits; 
titanium, nickel, cobalt and palladium silicides are common choices^i. We chose to study 
titanium silicide (TiSi 2 ) because of its technological importance in VLSI applicationsi^. Ti¬ 
tanium silicide (TiSi 2 ) exists in three phases: an orthorhombic base centered C49 phase, 
an orthorhombic face centered C54 phase and a hexagonal C40 phasel^^l^. TiSi 2 is exper¬ 
imentally synthesized by evaporating a thin layer of Ti on Si followed by annealing. The 
annealing temperature determines the structure of TiSi 2 phase; the C49 phase occurs when 
annealed at relatively low temperatures of about 700-870 K while the C54 phase forms at 
temperatures beyond 900 In the present work, a C49 titanium silicide (metal)- silicon 
(semiconductor) interface is chosen as a model metal-semiconductor interface to analyze 
thermal transport using hrst-principles DFT calculations. 

The remainder of the paper is organized as follows: we describe computational details 
in Section mi followed by a brief description of the mathematical formulation of electron- 
phonon coupling in Section IIIII In Section IIVI we present results from DFT calculations on 
bulk TiSi 2 and a TiSi 2 -Si interface. Complete spectral dispersion of electrons and phonons 
along with results for electron-phonon coupling are presented in Section HVl In Section |Vl we 
analyse the DFT results and estimate electron-phonon thermal conductances, and compare 
them with interfacial phonon conductance obtained from simplihed models. 
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II. COMPUTATIONAL METHODS 


All first-principles calculations in this paper are based on the DFT framework as imple¬ 
mented in Quantum Espresso^, which uses a planewave basis set and an ultrasoft pseudopo¬ 
tential to model the potential of nucleus and core electrons. The exchange correlation energy 
is approximated with a generalized gradient approximation (GGA) using the Perdew-Wang 
91 functional form, which has been used and well-tested in prior work on TiSijf^. An energy 
cutoff of 680 eV (50 Ry) is used to truncate the plane-wave expansion of the wavefunctions, 
and an energy cutoff of 6800 eV (500 Ry) is imposed on the plane-wave basis used to rep¬ 
resent charge density. The Brillouin zone is sampled with a uniform 10 x 10 x 10 grid of 
k-points in calculations (for the structural relaxation) of bulk TiSi 2 , and a grid of 10 x 10 x 1 
is used in the interface calculations. The cutoff energies and k-point grids are chosen to 
ensure that further increase in cutoff energies or refinement of k-point grids changed the 
total energy of the bulk and interface structures by less than 10 meV. 

A full relaxation (of both internal atomic positions and cell parameters) is performed for 
the bulk TiSi 2 calculations while a more restricted relaxation (see Section IIV Bp is performed 
for the TiSi 2 -Si supercell. The relaxation process terminates when the Hellmann-Feynman 
forces on each atom decrease below 0.025 eV/A, and when the change in energy between 
successive relaxation iterations is less than 1 meV. The band alignment and Schottky bar¬ 
rier height (SBH) are estimated using the electrostatic potential method and the projected 
density of states method. Both methods are found to give similar estimates of the SBH. 
The dynamical matrices (and phonons) are obtained using density functional perturbation 
theory (DFPT) for which the dynamical matrices are computed on a 2 x 2 x 2 grid for bulk 
G49 TiSi 2 and on a 2 x 2 x 1 grid for the supercell consisting of the TiSi 2 -Si interface. Us¬ 
ing this approach with Fourier interpolation, dynamical matrices are obtained at arbitrary 
q-vectors. 


III. MATHEMATICAL FORMULATION OF ELECTRON-PHONON 
COUPLING 

The definitions of electron-phonon coupling matrix elements, phonon linewidth and the 
associated Eliashberg function are reviewed below for the sake of completeness and clarity. 
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and are discussed in Refs. ll8l-l2Cll in more detail. The coupling matrix element for 

the scattering of an electron with wavevector k in the band z/ by a phonon of wavevector q 
and polarization p to a state with wavevector k + q in band fo is given by: 
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where 'ip denotes the electron wavefunction and 0 denotes the phonon eigenvector. The 
matrix element of electron-phonon coupling is given by the dot product between the phonon 
eigenvector and gradient of the self-consistent potential Vscf with respect to atomic displace¬ 
ments. Uqp is the phonon frequency, and M denotes the mass of atom (we hrst consider a 
case of mono-atomic unit cells). The formula for multi-atom unit cells involves a summa¬ 


tion over atoms in the unit cell, and is given in Ref. Il8|. The phonon linewidth (inverse of 


relaxation time) 7 qp due to electron-phonon scattering is dehned as: 
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where denotes the energy of electron in state fez/, Ef is the Fermi energy and VLbz is the 
volume of the Brillouin zone. The foregoing expression for phonon linewidth is valid for low 
temperatures at which electron-phonon interactions are conhned to a narrow region around 


the Fermi surface (see Ref. 


2ll ). The electron wavefunctions used to compute the coupling 


matrix elements are obtained on a 20 x 20 x 20 k-point grid in bulk TiSi 2 and on a 20 x 20 x 1 
k-point grid for the TiSi 2 -Si interface supercell. Such hne k-point grids are required for the 
convergence of the delta functions in Eq. (|2]). The Eliashberg function a^Eiu) is dehned in 
terms of a summation of the phonon linewidthsi^i^: 
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where D{Ef) is the electron density of states (DOS) at the Fermi energy. We now dehne a 
local Eliashberg function in a manner similar to the dehnition of the focal phonon DOS: 
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where I denotes atomic index within a nnit cell. a^Eiiu) is the local Eliashberg fnnction 
projected onto the atom I, and the snm over all the local (atomic) Eliashberg fnnctions gives 
the conventional total Eliashberg fnnction. 


TABLE I. Comparison of equilibrium lattice parameters of TiSi 2 in bulk C49 structure obtained 
from the present calculation with first-principles and experimental data reported in prior literature. 
The relaxed internal atomic positions (in the basis of the primitive lattice vectors) obtained from the 
present calculation are Ti: (±0.1,± 0.1,±0.25) and Si: (±0.75,±0.75,±0.25), (±0.44,±0.44,±0.25). 
VASP (LDA and GGA) values are obtained from Ref. Il7|. The range of experimental data is also 
obtained from Ref. [13 and other references cited therein. 


Method 

a (A) 

b(A) 

c(A) 

Present (QE, GGA) 

3.55 

13.62 

3.58 

VASP, LDA 

3.49 

13.27 

3.52 

VASP, GGA 

3.55 

13.62 

3.58 

Experiment 

3.55-3.62 

13.49-13.77 

3.55-3.65 


IV. ELECTRONIC STRUCTURE AND PHONONS OF BULK TiSis AND 
TiSia-Si INTERFACE 

A. Bulk C49 TiSia 


C49 TiSia exists in an orthorhombic structure with six atoms (2 Ti and 4 Si) per primitive 
unit cell, with primitive lattice vectors: 


^ a- b '• 

= 2“ + 2^ 
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FIG. 2. Electron bandstructure and DOS of bulk C49 TiSi 2 . 


Table U shows a comparison of the relaxed lattice parameters obtained from the present 
calculation with experimental and hrst-principles calculations reported in prior literature. 
The present results agree well with the GGA calculations of Ref. llT], and also fall within the 
typical DFT errors of experimental data. Our calculated electronic structure and density 
of states (see Figure |2]) agree very well with the results of linear augmented-plane-wave 
calculations of Mattheiss and HenseU^. 


Figure [3^ shows the phonon dispersion obtained from our DFPT calculations. We are not 
aware of any prior reports on the phonon dispersion of G49 TiSi 2 . However reports on the 
Raman characterization of TiSi 2 are available in the literature^*^. We hnd excellent agree¬ 
ment (about 1% error) of the DFPT results with the experimental Raman modes. These 
results for lattice structure and bandstructures of bulk G49 TiSi 2 validate the accuracy of 
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FIG. 3. a) Phonon dispersion of bulk C49 TiSi 2 . The squares correspond to Raman-active modes 
predicted from DFPT calculations, and the horizontal lines correspond to experimental Raman 
frequency shifts characteristic of the C49 phase^^. b) Phonon DOS and Eliashberg function of bulk 
C49 TiSis. 


calculational parameters, the exchange correlation energy functionals and pseudopotentials 
used in the present DFT formulation. 

The range of the Eliashberg function of TiSi 2 (see Figure (Sb) is of the same order as that 
for typical metals such as Al, An and Nai^. The peaks of the Eliashberg function closely 
follow the peaks of the phonon DOS, as expected from the form of Eq. (I3]). Intuitively, the 
Eliashberg function can be visualized as the product of phonon DOS and electron-phonon 
coupling linewidths. 
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FIG. 4. Two interface configurations between (0 0 1) Si and (0 1 0) C49 TiSi 2 . Blue spheres are Si 
atoms and magenta spheres are Ti atoms, a) Interface Configuration 1 with 28 atoms in the unit 
cell. The red dotted boxed region are the atomic layers considered in the calculation of Eliashberg 
function for interface modes in Section |Vl b) Interface Configuration 2 with 32 atoms in the unit 
cell. 

B. C49 TiSi 2 -Si Interface 
1. Structure and Stability 

A number of experimental studies on the epitaxial growth of C49 TiSi 2 on Si have been 
reported in the literature^^^— , however with a wide range of the Miller indices of TiSi 2 and 
Si planes that form the interface. Multiple possible epitaxial orientations have been reported 
within the same sample based on detailed TEM and diffraction studies of the interface. Wang 
et al.-^ performed a detailed study on the extent of in-plane lattice mismatch for different 
Miller indices and relative orientations of the TiSi 2 and Si planes forming the interface. 
Among these, the epitaxy of the (0 1 0) TiSi 2 plane parallel to the (0 0 1) Si plane with 
the [1 0 0] direction in TiSi 2 parallel to the [1 1 0] direction in Si had the smallest lattice 
mismatch in terms of superlattice area. Hence, this set of planes at the interface are likely to 
be the most probable in terms of energy cost of lattice mismatch and has also been reported 
to be the most prominent epitaxy observed^^*^^. All further studies in this paper have been 
performed for this particular orientation at the interface. 

In all DFT simulations of the interface, TiSi 2 is stretched along the [1 0 0] and [0 0 1] 
directions to match the bulk lattice constant of silicon. This condition mimics the experi¬ 
mental synthesis process where a thin him of TiSi 2 is typically grown on a bulk or thick Si 
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substrate. Hence, Si can be assumed to maintain its bulk in-plane lattice constant while the 
thin him of TiSi 2 stretches in the plane of the interface. Si has a bulk lattice constant of 3.84 
A cubic/ in the (0 0 1) plane while TiSi 2 has bulk lattice constants of 3.55 A and 3.57 
A (see Table H]) along the (0 1 0) plane. Hence a tensile strain of 7% is required in the (0 1 
0) plane of TiSi 2 to match the bulk lattice constant of Si. Even though this is a signihcant 
strain, we chose to work with this interface for the reasons explained above. Also, this par¬ 
ticular epitaxial structure exhibits the smallest supercell in-plane area and hence reduces the 
number of atoms needed to form the supercell. This feature is important for computational 
efficiency in the subsequent phonon and electron-phonon coupling calculations. 

Figure 0] shows unit cells of two interfacial conhgurations studied in the present work. 
Each supercell has two Si-TiSi 2 interfaces: one near the center and the other near the 
right edge of the supercell. The supercells are constructed such that both interfaces are 
structurally identical in the unrelaxed conhguration. The atomic structures in Figure 0] are 
relaxed by keeping the in-plane lattice constants hxed to the bulk lattice constant of Si. 
All internal atomic positions and the length of the unit cell along the transport (out-of¬ 
plane) direction are relaxed. A contraction in the length of the unit cell along the transport 
direction is observed after relaxation to minimum energy for both interfacial conhgurations, 
as expected from the tensile strain imposed on TiSi 2 along the in-plane direction. 

The interface energies are determined (see Table [TTl) using the following formula: 

rp _ rps _ 

„ _ ^supercell ^Si ^TiSi2 ( 

-Winter face 

where Eguperceii is the energy of the supercell shown in Figure 0] after relaxation. Ffj is calcu¬ 
lated by relaxation after removing all atoms belonging to TiSi 2 in the supercell, and 
is calculated after removing all atoms belonging to Si (the superscript s indicates a surface). 
The factor of two accounts for the presence of two interfaces within a supercell. Both super¬ 
cells in Figure 0] are constructed such that the Si plane of TiSi 2 forms the interface with Si. 
TiSi 2 in the C49 structure has two kinds of Si atomsl^, and we obtain two conhgurations of 
the interface at the planes of these two types of Si atoms. The area normalized interfacial 
energy is given as Einterface/^ where A is the cross-sectional area. Both interfacial conhgura- 
tions in Figure0]have similar interface energies of about 0.2 eV/A . Table UTl also shows the 
normal stresses along the Cartesian directions after relaxation. High negative stresses along 
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the in-plane directions {x and y) are expected due to the tensile strain imposed by stretching 
the TiSi 2 lattice. The stress along the transport direction vanishes upon relaxation of the 
length of the supercell in the z-direction. 


TABLE II. Properties of the two interface configurations shown in Figure 01 


Conhguration 

Interface Energy (eV/A^) 

Stress (kbar) 

p-type SBH (eV) 

1 

-0.18 

o-xx = -54.2, ayy = -70.23, a^z = -0.12 

0.47 

2 

-0.19 

o'xx = —65.68, ayy = —83.35, a^z = —0.06 

-0.29 


2. Band Alignment 


Table HIl also shows the p-type SBH of the two interface conhgurations calculated using 


the band line-up method of Ref. 


271 . According to this method, the p-type Schottky barrier 


Esbh is calculated using the following formula: 


Esbh = (Bj'®" - - KL) + AC (7) 

where = 7.62 eV is the difference between the Fermi energy and the macro- 

scopically averaged electrostatic potential in bulk TiSi 2 . is computed for 

a strained bulk TiSi 2 crystal with the same in-plane strain that is used to construct the 
lattice-matched interface. The Fermi energy in Figure |5] is located at an energy level of 
7.62 eV above the macroscopically averaged electrostatic potential on the TiSi 2 side of the 
interface. The foregoing method of locating the Fermi level gives almost the same result 
as that obtained directly from the analysis of the self-consistent DFT solution of the re¬ 
laxed supercell. Ey^j^ — = 4.81 eV is the valence band maximum (VBM) calculated 

with reference to the macroscopically averaged electrostatic potential in the bulk Si crystal. 
Although DFT is known to under-predict the band gap, the energy of VBM and its offset 
relative to Fermi energy are expected to be quite accurate. Hence this value is used to locate 
the VBM on the Si side of the interface in Figure [5l The conduction band minimum (CBM) 
is located using the experimental^^ band gap (1.1 eV) of bulk Si. AV denotes the height of 
the step in the total potential at the interface and is the difference between macroscopically 
averaged electrostatic potentials on either side of the interface (see Figure |5]). In interface 
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FIG. 5. Planar average electrostatic potential (blue continuous) and macroscopic average electro¬ 
static potential (red dashed) plotted along the transport direction of interface Conhgurations 1 
and 2 shown in Figure 01 The VBM in Si and Fermi energy in TiSi 2 are located with reference to 
the respective macroscopic electrostatic potentials using separate calculations on bulk Si and bulk 
TiSia. 


Configuration 1, the Fermi energy is found to lie between the valence and conduction bands 
of Si, and a p-type Schottky barrier of 0.47 eV is predicted. Our estimate of 0.47 eV is 
consistent with experimental SBH measurements of 0.5 - 0.6 eV in Ref. 
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However in the 


interface Conhguration 2, the Fermi energy lies 0.29 eV below the VBM of Si, and hence it 
behaves as an ohmic contact. 

The estimates of SBH using the electrostatic potential method are also conhrmed using 
the projected density of states (PDOS) method. Figure [6] shows the PDOS on a Si atom 
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FIG. 6. PDOS on a Si atom farthest from the interfaces in the supercell. Configurations 1 and 2 
refer to the two interface structures shown in Figured! 

that is farthest from the interfaces in the supercell. This PDOS is expected to mimic the 
bulk DOS of Si. The VBM is located about 0.6 eV below the Fermi energy of interface 
Conhguration 1 while the VBM is located about 0.25 eV above the Fermi energy of interface 
Conhguration 2. These values are reasonably close to the SBH estimates obtained from the 
electrostatic potential line-up method. The potential line-up estimates are however more 
reliable because the spatial convergence of local DOS to the bulk value is slower than the 
spatial convergence of potentials or charge density^l. The foregoing results illustrate that 
the details of the atomic scale structure of an interface decide the behavior of a metal- 
semiconductor contact. 


3. Phonons and Electron-Phonon Coupling 

DFPT calculations are computationally intensive, and hence we focus only on interface 
Conhguration 1 (Schottky interface) for further phonon and electron-phonon coupling cal¬ 
culations. Figure [7| shows the phonon dispersion of the supercell with TiSi 2 -Si interface 
Conhguration 1 calculated in the two-dimensional Brillouin zone. Real vibrational frequen¬ 
cies suggest that the interface conhguration is stable in spite of the large tensile strain 
imposed on TiSi 2 to match its lattice with Si. We also observe the existence of many soft 
(low frequency and dispersionless) interfacial modes that are expected to further stabilize 
the structure at higher temperatures. 
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FIG. 7. Phonon dispersion of interface supercell Configuration 1. 


We now examine the phonon DOS of TiSi 2 -Si interface and compare it with that of bulk 
TiSi 2 and bulk Si (Figure [H]). The phonon DOS for bulk TiSi 2 and Si have been scaled 
by factors of 2 and 8 respectively because the interface supercell (see Figure IHi) contains 
2 unit cells of TiSi 2 (4 Ti and 8 Si atoms) and 8 unit cells of Si (16 Si atoms). The bulk 
phonon DOS of TiSi 2 is plotted both for the relaxed unstrained crystal and for a structure 
with 7% tensile strain. We observe that the peak DOS in bulk TiSi 2 (both strained and 
unstrained) occurs at frequencies in the range of 200-300 cm“^ while that of bulk Si occurs in 
the frequency range of 150-200 cm“^. The interfacial phonon DOS has a relatively broader 
peak between 150-300 cm“^ that encompasses both these frequency ranges. This result 
suggests that the interfacial phonon DOS follows an average of the bulk structures forming 
the interface. Similar shifts in the DOS of interfacial modes have been reported for In/Si 


interfaces in Ref. 


30 


A comparison between the Eliashberg function (Figure [9]) of the supercell with interface 
Conhguration 1 with that of strained and unstrained bulk TiSi 2 indicates that the average 
Eliashberg function for the interfacial structure is of the same order of magnitude as that 
for bulk TiSi 2 . Since TiSi 2 in the interfacial structure is strained by 7%, Eliashberg function 
of the interfacial configuration follows the Eliashberg function of strained bulk TiSi 2 . The 
tensile strain in the structure also produces a distinct red-shift in the low-frequency peaks of 
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FIG. 8. Phonon DOS of the TiSi 2 -Si interface supercell Configuration 1 compared with that of 
bulk TiSi 2 and bulk Si. 



FIG. 9. Eliashberg function of the TiSi 2 -Si interface Configuration 1 compared with that of bulk 
TiSi2. 


the Eliashberg function. The Eliashberg function of the strained crystal peaks in the range 
of 150-200 cm“^ while that of unstrained bulk TiSi 2 peaks around 300 cm“^. 

The Eliashberg function can be decomposed spatially using the phonon eigenvectors as 
shown in Eq. (jl]). Figure flOh shows the spatially decomposed Eliashberg function for three 
different regions of the supercell: atoms belonging to Si that are farthest from the interface, 
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atoms belonging to TiSi 2 that are farthest from the interface and atoms near the interfacial 
region. The Eliashberg function is summed over 6 atomic layers in each of these three regions. 
As expected, the strength of electron-phonon coupling is strongest on the metal side of the 
interface and weakest on the semiconductor side. The interfacial contribution to Eliashberg 
function is dominant in the 150-200 cm“^ range of frequencies (the modes red-shifted due to 
the tensile strain). Figure ITOb shows the Eliashberg function projected on different atomic 
planes along the supercell for two specihc phonons with frequencies of 170 and 250 cm“^. 
The Eliashberg function is weak on the Si side but increases sharply across the interface 
to approximately ten times its value on the Si side. The decay in the Eliashberg function 
beyond layer 25 is due to the second interface (and the Si region beyond it) that exists 
because of the translational periodicity used in the calculations. The Eliashberg function 
on the Si side of the interface suggests that the phonon mode at 170 cm“^ couples more 
strongly with the electronic states penetrating into the semiconductor than the mode at 250 
cm“^. As a result, this mode is expected to contribute more actively to the cross-interface 
electron-phonon thermal conductance. Also, the local Eliashberg function on the silicon side 
of the interface does not decay immediately to zero as the local electron DOS at the Fermi 
energy is also not exactly zero for a silicon atom farthest from the interface and the edges of 
the supercell (see PDOS in Figure [6]). This is due to the slow convergence of local electron 
DOS to the bulk Si DOS at such small distances (~10 A) from the interface layer. 


V. THERMAL CONDUCTANCE OF THE TiSis-Si INTERFACE 

In this section, we report predictions of thermal conductance along the different heat 
transfer pathways (electron-phonon, phonon-phonon) shown schematically in Figure [H The 
objective is to obtain a quantitative estimate of the contributions of each of these pathway^ 
which could help assess the validity of simplihed models such as those proposed in Refs, [sl 
andO. We consider three different contributions to thermal conductance (see Figure [T]) in 
this section: 


• Gep,b is the electron-phonon thermal conductance in the bulk region of the metal and 
is obtained from the Eliashberg function of strained bulk TiSi 2 . This conductance is 
used in the models of Refs. ^ andO (Mechanism A in Figured]). 


17 





Distance along supercell (Angstrom) 

FIG. 10. a) Spectral variation of the Eliashberg function over three different regions of the supercell 
with interface Configuration 1. b) Variation of the Eliashberg function along the transport direction 
of the interface Configuration 1 for phonon frequencies of 170 and 250 cm~^. 

• Gep^i is the electron-phonon thermal conductance contributed by degrees of freedom 
near the interface region and physically represents the coupling of electrons in metal 
with the interfacial or joint phonon modes. Gep^i is obtained from the Eliashberg 
function of the interface supercell. This conductance was proposed to be important in 
the thermal conductance of a Pb/diamond interface^ and was also considered in the 
works of Sergeev^i^. Since the joint phonon modes exist on both sides of the interface, a 
part of the energy transferred from electrons to interfacial phonon modes is transferred 
across the interface. Hence, this conductance relates to direct coupling of electrons to 
phonons across the interface but is often neglected with the assumption that electrons 
in metal are thermally insulated from the semiconductor at the interface^'^'^. 
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• Gpp is the phonon thermal conductance of the interface and can be evaluated using 
advanced techniques such as the atomistic Green’s function method^^. In this work, 
we adopt a simplihed approach to estimate Gpp using the diffuse mismatch model^ and 
also obtain an upper bound from the radiation limit. 


The energy transfer rate (per unit volume) J^p from electrons to phonons can be expressed 
in terms of the Eliashberg function using the following expression (see Ref. 2 for a detailed 
derivation): 


Jep = 27rD{Ef) j {huj)WF{u){rsEi^,Te) - rp,j,{u,T,))du ( 8 ) 

0 

where and T; denote the electron and lattice temperatures respectively, and de¬ 
notes the Bose-Einstein distribution function. The electron-phonon conductance G^p can be 
dehned for a small difference between the electron and lattice temperatures as: 

OO 

G,p = 2nD{Ef) j {hufa^E{u)^^du (9) 

0 


This expression is used to compute the electron-phonon conductances Gep,b and Gep,i where 
the hrst-principles Eliashberg function of bulk TiSi 2 and the interface supercell are used 
in the integrand. The latter (Ggp,*) is computed by summing the local Eliashberg function 
projected on 19 atomic planes around the interface near the center of the supercell. The 
contribution to the Eliashberg function from other layers near the edges of the supercell are 
not considered in order to extract only the contribution of electron-phonon coupling near 
a single interface. These 19 atomic layers are shown in Figure 0^ and correspond to the 
atomic index I (see Eq. (|3|)) in the range of 5 to 23 where Z = 15 is the location of the 
interface. The electron-phonon conductance computed from Eq. ([9]) has units of W/m^K 
and needs to be converted to W/m^K to obtain an area-normalized conductance that can 
be compared with Gpp. We do this by multiplying by a typical length scale for electron- 
phonon energy transfer. A length scale of 3Aep which corresponds to 95% energy loss was 
suggested in Ref. ^ where Xep is the electron mean free path due to scattering from phonons 

& determined electrical resistivity (p) and carrier 
, Xep can be estimated using the simple formula 


in bulk TiSi 2 . Using the experimenta’ 
concentration (n) of C49 TiSi 2 in Ref. 
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\ep = mvf/ne^p. The low-temperature residual resistivity is subtracted from the actual 
resistivity at room temperature to obtain the resistivity due to electron-phonon scattering 
alone, excluding other scattering effects such as grain boundary and impurity scattering. 
The Fermi velocity of bulk C49 TiSi 2 is Vf = 0.49 x 10® rn/sl^. We obtain an electron mean 
free path due to phonon scattering of about 4 nm at room temperature and use a length scale 
of 10 nm for converting the volumetric electron-phonon conductance to an area-normalized 
conductance. 


The spatial extent of the joint phonon modes depends on the anharmonicity of the local 
interfacial structure and is computationally expensive to predict from hrst-principles calcu¬ 
lations. A heuristic approximation for the spatial extent of joint/interfacial phonon modes 
was provided by Huberman and Overhauser- who assumed that the joint modes extend to 
a distance of the bulk phonon mean free path on either side of the interface. We estimate 
the phonon mean free path Ap in bulk TiSi 2 using the kinetic theory expression for ther- 


le phonon contribution to thermal conductivity Kp 
by subtracting the electronic thermal conductivity 


mal conductivity Kp = {l/3)cyVgAp. T' 
was estimated as 3.4 W/mK in Ref. l34 
(obtained using the Wiedemann Franz law) from the total thermal conductivity measured 
experimentally. An average phonon group velocity of 5900 m/s is estimated using the phonon 
dispersion in Figure The specihc heat capacity of C54 TiSi 2 was experimentally mea¬ 


sured and reported as 66.9 J/mol K in Ref. 
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We assume the same value for C49 TiSi 2 as 


the specihc heat is not expected to be signihcantly different. This gives a phonon mean free 
path in C49 TiSi 2 as 7 A while the phonon mean free path in silicon is about 400 A—. Such 
a large difference in phonon mean free paths indicates that the energy of the joint modes is 
predominantly concentrated on the Si side of the interface. Hence the electron-joint phonon 
mode conductance Gep.i calculated using the Eliashberg function of the supercell containing 
the interface is expected to be a good estimate of the cross-interface electron-phonon thermal 
conductance. We note that this approach is similar to the Pb/diamond interface considered 


in Ref. 


a 


where the energy of the joint modes is concentrated primarily in diamond due to 
the small phonon mean free path in Pb. 

The phonon-phonon interfacial conductance Gpp is estimated from the Landauer formula: 


G 


pp 



hu{q,p)vg^^ 


^Ibe 

dT 


T{uj{q,p))£q 


( 10 ) 


q,gz>0 
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FIG. 11. Thermal conductance due to electron-phonon coupling for bulk TiSi 2 (Gep^b) and the 
interface (Gep^i) along with interfacial phonon-phonon thermal conductances from the radiation 
limit and DMM. The arrows indicate the y-axis to which each curve corresponds to. 


where Vg^z denotes the phonon gronp velocity along the transport direction. The trans¬ 
mission fnnction T{oj{q,p)) is estimated in the harmonic radiation limit (npper limit for 
harmonic interface conductance) and also from the diffuse mismatch model (DMM) using 


the calculated full phonon dispersions of strained bulk TiSi 2 anc 


details on the radiation limit of interface conductance, and Ref. 
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bulk Si. Ref. y provides 
describes implementation 


of DMM using full phonon dispersion rather than the more common Debye approximation. 

Figure [TT] shows the variation of Gpp, Gep,b and G^p^i with temperature. While Gep,i is 
lower than Gpp and Gep,b at room temperature it is not low enough that direct coupling from 


electrons in metal to interfacial p 


series such as that used in Refs. 


lonon modes can be neglected. A model of resistors in 


and ^ assumes electron-phonon coupling across the 
interface to be insignihcant compared to electron-phonon coupling in the bulk of the metal. 
Our results however indicate that the coupling strength of electrons to the interface phonon 
modes at a TiSi 2 -Si interface is of the same order of magnitude as the electron-phonon 
coupling in bulk TiSi 2 metal. Thus, a simple series resistor model may not be a valid 
approximation in general to the thermal interface conductance. Also, the electron-phonon 
thermal conductances (both G^p^ and G^p^i) are of the same order of magnitude as the 
phonon-phonon conductance across the TiSi 2 -Si interface at room temperature. Moreover, 
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we note that the electron-phonon interaction across an interface provides a direct coupling 
of the primary thermal carriers in the respective materials. Therefore, this coupling process 
should be considered a parallel path to the all-phonon cross-interface transport process. 

Finally, we note that although the results on electron-phonon coupling have been pre¬ 
sented only for interface Conhguration 1, the strength of cross-interface electron-phonon 
coupling is expected to be even stronger for interface Conhguration 2 because the latter 
forms an ohmic contact and hence electrons from the metal experience no energy barrier for 
coupling with phonons in the semiconductor. These hndings contrast with those of Ref. 
which concluded that electron-phonon coupling resistance can be ignored at temperatures 
above 200 K, albeit for different metal-semiconductor interfaces than considered in this work. 


VI. CONCLUSIONS 

The contribution of electron-phonon coupling, especially the direct coupling of electrons 
in a metal to phonons in a semiconductor, is often neglected in the study of thermal transport 
across a metal-semiconductor interface. The present ab initio studies on a technologically 
important TiSi 2 (metal) - Si (semiconductor) interface however indicate that the strength 
of coupling between electrons and interfacial joint phonon modes is of the same order of 
magnitude as the strength of electron-phonon coupling in bulk metal. Further, the thermal 
conductance from the coupling of electrons to interfacial or bulk metal phonons is compa¬ 
rable to the phonon-phonon conductance across the interface. These results highlight the 
importance of considering the direct coupling of electrons in metal to phonons in semiconduc¬ 
tor as an alternate pathway for interfacial thermal transport in addition to the conventional 
phonon (metal) - phonon (semiconductor) energy transfer pathway. 

Looking toward future developments, we hrst note that the non-equilibrium Green’s func¬ 
tion (NEGF) method for electron transport and the analogous atomistic Green’s function 
method for phonon transport are well established in contemporary literature. However, re¬ 
ports of coupled electron-phonon transport using the Green’s function method are scarce 
and exist only for simple systems such as ID atom chains^i^. We are currently working on 
the development of a NEGF-based electron-phonon coupled transport solver to which the re¬ 
sults from hrst principles calculations (such as inter-atomic force constants, electron-phonon 
coupling matrix elements and the associated Eliashberg function) are direct inputs. 
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